In-medium meson effects on the equation of state of hot and dense nuclear matter 
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The influence of the in-medium mesons on the efi'ective nucleon mass and in turn on the equation 
of state of hot/dense nuclear matter is discussed in the Walecka modeL Due to the self-consistent 
treatment of couplings between nucleons and a and u mesons, the temperature and density de- 
pendence of the efi'ective hadron masses approaches more towards the Brown-Rho scaling law, and 
the compression modulus K is reduced from 550 MeV in mean field theory to an accepted value 
318.2 MeV. 
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The study of strongly interacting nuclear matter un- 
der extreme condition realized in relativistic heavy ion 
collisions has attracted a lot of attention during recent 
years(for example, see Ref. [1] and references therein). In 
high energy nuclear collisions the temperature may ap- 
proach to several hundred MeV's or the density to sev- 
eral times of the normal nuclear density pQ. It is believed 
that the property of hadrons in such high temperature 
or density will be quite different from that in vacuum. 
For instance, the behaviors of the in-medium hadronic 
masses may be qualitatively expressed by the so-called 
Brown-Rho scaling law [2]. 

The quantum hadrodynamics QHD-I model (Walecka 
model ) and its extensions have been widely used to dis- 
cuss the property of the symmetric nuclear matter and 
finite nuclei [3-7]. For example, the saturation proper- 
ties at the normal density can be explained successfully 
by the mean field approximation. 

In the Walecka model, the nucleons interact with each 
other through the exchange of a and lu mesons. The a 
exchange gives the attractive force, while the lo exchange 
the repulsive force. The mean field theory (MFT) has 
been widely used to discuss the nuclear matter satura- 
tion property at the normal density, while it was found 
that the obtained compression modulus K ~ 550 MeV is 
larger than the acceptable value 300 ± 50 MeV [8] . The 
reason for so high compression modulus is that the effec- 
tive nucleon mass drops too fast with the density. By tak- 
ing into account the vacuum fluctuations, the behavior 
of the effective nucleon mass with density may be cured 
to some extent, but the obtained compression modulus 
K ~ 450 MeV is still larger than the acceptable value. 
To solve this problem, the a — lu model including non- 
linear interaction terms ba^ + ca^ of a and the Zimanyi- 
Moszkowski (ZM) model have been proposed [6,7]. With 
more adjustable parameters in the nonlinear a — to ver- 
sion, the effective nucleon mass drops slowly and the com- 
pression modulus K may be in the region of "appropriate 
values." The behaviors of the effective nucleon mass and 



the compression modulus K may be also cured to some 
extent with the nonrenormalizable ZM model [7]. 

On the other hand, the in-medium effects of hadrons 
under extreme environment have been emphasized in the 
QHD-like framework. Since Walecka model's Lagrangian 
is a truncated form of a chirally symmetric Lagrangian 
with the vector and scalar fields taken as chiral-singlet 
fields [9,10], it is argued that the predictions about the 
effective masses of hadrons in hot / dense environments by 
QHD are intrinsically consistent with the chiral symme- 
try restoration [11,12]. This simple model includes im- 
plicitly vacuum effects and subnucleon structures. With 
QHD-I and by considering the Dirac sea contribution, 
the obtained effective masses of light vector mesons drop 
down with increasing density/temperature. However, in 
discussing the hadronic masses, the relevant equations of 
QHD-I model had not been solved with a consistent man- 
ner, i.e., first solving the equations for the effective nu- 
cleon mass and chemical potential fi*j^ (for finite tem- 
perature occasion) in MFT or RHA (relativistic Hartree 
approximation) approach and then using the obtained 
and /x^ to get the full propagators for mesons a and 
uj and even p [13-16]. Therefore, the in-medium meson 
effects were not reflected in the effective nucleon mass as 
pointed out by Bhattacharyya et al. in Ref. [17]. Our 
motivation is to study the effect of nonperturbative in- 
medium modification of a and lo on the effective nucleon 
mass and on the equation of state (EOS) of hot/dense 
nuclear matter by using renormalizable original version 
QHD-I. With Dyson-Schwinger Green function approach, 
the in-medium resummed nucleon and meson propaga- 
tors are treated self-consistently. 

The full description for the Lagrangian of Walecka 
model can be found in Refs. [4,18]. Under the mean 
field approximation, the full nucleon propagator in the 
medium as indicated in Fig.l is attributed to the calcu- 
lation of the tadpole self-energies with finite temperature 
field theory: 
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where po = (2n + l)7rTi + with T and being tem- 
perature and effective baryon chemical potential, respec- 
tively, and the symbol /p = / (^J^s ■ With the residue 
theorem, one can divide the self-energies into the vac- 
uum fluctuation and the matter (obviously related to the 
distribution functions) parts, and define the effective nu- 
cleon mass and effective chemical potential /x^ [4] , 
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In Eqs. (3), (6) and (7), the explicit expressions for the 
vacuum fluctuation contributions AM*^^ and Ae*^^ can 
be found in Refs. [4,19]. The thermodynamics compen- 
satory term B* in Eq. (7) can be determined uniquely 
by the thermodynamics self-consistency relation between 
the energy density e and pressure p. For MFT and RHA 
approaches, this term is vanishing. 

The properties of mesons are studied by the meson 
propagators in medium, which are normally calculated by 
using the random phase approximation (RPA) with the 
full nucleon propagator in Fig.l. The Dyson-Schwinger 
equation of the uo meson propagator is indicated in Fig. 2 
and its solution for the full propagator D^^" can be de- 
termined by the polarization tensor 11'*'^ (fc). 
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where the baryon density 
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The spin-isospin degenerate factor is 7 = 4 for symmetric 
nuclear matter and 2 for neutron star. 
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FIG. 1. Diagrammatic representation for the nucleon prop- 
agator in relativistic Hartree approximation (RHA). 



where is the bare propagator. Using the Feynman 
rules at finite temperature [20], one has 



W^{k) = 2glTY, f Tr 



N 



• (9) 
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FIG. 2. The in-medium uj propagator in random phase ap- 
proximation. 

Different from the occasion in vacuum, there are two 
independent elements 11^ '-^^ for the polarization tensor 

n^''(fc) in medium, 

u'^'ik) = n^(fc)pr + n^(A;)Pr> 
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One can see that the self-consistent equation of effec- 
tive nucleon mass in RHA is just the result in the MFT 
plus the vacuum fluctuation contribution AM*^^. Simi- 
larly, the energy density and the pressure are the MFT 
results plus the vacuum fluctuation contribution Ae*^^, 
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where the Pj^^ and Pj,^ are the standard projection ten- 
sors [20]. The ingredients JlL{k) and Tirik) are deter- 
mined by 

iiL{k) = ^noo(fc), nT(fc) = \p^Tii^{k). (11) 

The pole position of the full propagator D^^" determines 
completely the dispersion relation of u) meson excitations 
in medium. The expressions for the various components 
of n'^''(fc) are similar to those of p as given in Refs. [18,21] 
except the vanishing of the tensor coupling constant for 
Lo meson, i.e., = 0. 

Analogously to the discussion of uj meson, one can dis- 
cuss the property of a in the medium. Its diagrammatic 
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Dyson-Schwinger equation is similar to Fig. 
self-energy is [21,22] 
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There arc two kinds of effective meson masses related 
to our work. One is the pole mass m*, which can be 
matched to the results obtained by such as QCD sum 
rules. It is defined by the pole position of the meson 
propagator in medium by taking the limit |k| — > of 

n(fc), 
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The other is the off-shell mass m in the medium deter- 
mined by 



= + lim lim Il(k). 

fco^o|kHo 



(14) 



Due to the fact that the longitudinal and transverse com- 
ponents of the to polarization tensor 11^'' (fc) coincide with 
each other in taking the limit |k| — »■ 0, the subscripts L 
and T for Ilu;{k) are omitted. 

As calculated in the following, the meson property will 
be quite different from the vacuum scenario. To consider 
the back interactions of a and lu mesons with nucleon N 
from the point of view of self-consistency, the in-medium 
a and u propagators with vanishing four-momentum 
transfer should be used in determining the nucleon prop- 
agator as indicated by Fig. 1 due to the adapted Hartree 
approximation. Therefore, the meson masses rricr and m,^ 
in the effective nucleon mass and chemical potential equa- 
tions (3) and (4), and in the EOS's (6) and (7) should 
be the off-shell meson masses to's determined by Eqs. 
(14). It means that Eqs. (3) and (4) obtained by the 
RHA of the nucleon propagator and Eqs. (14) obtained 
by the RPA of the meson propagator form a closed set of 
equations, and should be solved simultaneously. 
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FIG. 3. (a) The binding energy for nuclear matter as a 
function of Fermi momentum at T = 0. The dot-dashed, 
dashed, and solid lines represent the MFT, RHA, and our re- 
sults, respectively; (b) The pressure as a function of rescaled 
density ps /po for two temperatures T = 20 MeV and 30 MeV. 

Considering the in-medium effects of the off-shell 

masses (14) on the EOS, one must at first refit the bind- 
ing energy at the normal nuclear density at zero tempera- 
ture by adjusting the coupling constants goj{a) as done in 
MFT and RHA. The coupling constants are listed in Ta- 
ble I, which are smaller than those of MFT and RHA. It 
is interesting to note that the EOS becomes much softer 
and the compression modulus 
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is reduced to = 318.2 MeV, which is acceptable in 
dealing with realistic nuclear matter. In Eq. (15), ejv 
is the binding energy {e/ps — Mj^). Of course, this in- 
teresting numerical result can be also obtained with the 
nonlinear a — lu version by adjusting the additional pa- 
rameters due to the renormalizability of original QHD-I 
as pointed out in the introduction [5]. However, it is 
interesting enough that it has been obtained here from 
the point of view of in-medium meson contribution on 
the nucleon and bulk property of nuclear matter. To the 
best of our knowledge, a similar numerical result of com- 
pression modulus K has also been obtained long ago by 
Ji in Ref. [23] with QHD-I through taking into account 
the polarization effects at zero temperature. It should be 
pointed out that our numerical results (especially ,9^; 
and K) are also different from those of Ref. [17] for the 
zero-temperature scenario. 
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One can further study the energy density and pressure. 
One should note that the consideration of the back in- 
teractions of in-medium mesons with nucleons will lead 
to the nonvanishing compensatory term B* in Eq. (7) 
due to the implicit higher order contribution from the 
resummed meson propagators instead of frozen ones. In 
principle, this term can be expanded into a series of 
— Mjv. With the energy density Eq. (6) and thermo- 
dynamics self-consistency condition at zero temperature 
p = p% {d/dps) {e/pb), this compensatory term can be 
fixed uniquely. In the upper panel (a) of Fig. 3, we give 
out the binding energy at T = 0, from which one can 
confirm that the equation of state really becomes softer 
and the compression modulus is much smaller than the 
results of MET and RHA. The first order liquid-gas phase 
transition in the low temperature case still exists as in- 
dicated by Eig.3(b) and the critical temperature Tc is 
about 21 MeV. 
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FIG. 4. The effective masses (solid lines ), (dotted 
lines) and m* (dot-dashed lines) as functions of rescaled den- 
sity pb/p» at T = 100 MeV (a) and temperature T at ps = 
(b), respectively. 

The pole masses Mt^, m* and m* determined by Eqs. 
(3) and (13) arc indicated in Fig. 4. Because of con- 
sidering the back interactions of a and uj mesons with 
nucleons, their temperature and density dependence ap- 
proaches more towards the Brown-Rho scaling law com- 
pared with the results of MFT/RHA. The key point is 
that the exchanged mesons are not bare but in-medium 
ones, their masses are self-consistently determined by the 
Dyson- Schwinger equations. The temperature and den- 
sity dependence of hadron masses studied here is analo- 
gous to that of hidden local symmetry theory of Harada 
and Yamawaki where such dependence is required by the 
Wilsonian matching to QCD [24]. It is interesting to 
note that our results are consistent with those obtained 
in recent works [25,26], where the Brown-Rho scaling fits 
naturally into the relevant framework. This consistency 
is attributed to the hidden chiral symmetry in QHD-I 
and should contribute to the understanding of hadronic 
matter under extreme conditions within the relativistic 
nuclear theory framework. 

In summary, we have investigated the in-medium me- 
son contribution to the EOS of hot/dense nuclear matter 
with the original Walecka model by treating the coupled 
nucleon and meson propagators self-consistently. This 
kind of contribution makes the EOS softer than the RHA 
and MET approaches. Due to the couplings between nu- 
cleon and in-medium mesons, the compression modulus 
K drops from about 550 MeV to about 320 MeV and 
the effective hadronic masses approach more towards the 
Brown-Rho scaling law, and the low temperature liquid- 
gas phase transition still exists. 
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TABLE I. The coupling constants ga and gu,, the effective nucleon mass and the compression modulus K at normal 
nuclear density po = 0.1484fm~* at T = 0. In MFT, RHA, and our approach(labeled as RHA+RPA), the nucleon and oj meson 

masses in vacuum arc taken to be A4n = 939 MeV, = 783 MeV, and the a meson mass in the vacuum is selected to be 
m<T = 520 MeV (in the literature with MFT) and 458 MeV. To compare with the results in MFT and RHA, we have also shown 
the medium dependent coupling constants Cg = g^{M% /fhi) and = gt{M% /fhtj). 
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MFT 


109.94 


191.04 


520 


358.49 


274.76 


547.2 


0.540 




85.286 


191.04 


458 


358.49 


274.76 


547.2 


0.540 


RHA 


54.289 


102.76 


458 


228.19 


147.78 


452.6 


0.731 


RHA-I-RPA 


42.229 


69.729 


458 


152.59 


86.263 


318.2 


0.803 
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